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Precise optical properties of metals are very important for accurate prediction of the Casimir force 
acting between two metallic plates. Therefore we measured ellipsometrically the optical responses 
of Au films in a wide range of wavelengths from 0.14 pan to 33 [im. The films at various thickness 
were deposited at different conditions on silicon or mica substrates. Considerable variation of the 
frequency dependent dielectric function from sample to sample was found. Detailed analysis of 
the dielectric functions was performed to check the Kramers-Kronig consistency, and extract the 
Drude parameters of the films. It was found that the plasma frequency varies in the range from 
6.8 eV to 8.4 eV. It is suggested that this variation is related with the film density. X-ray reflectivity 
measurements support qualitatively this conclusion. The Casimir force is evaluated for the dielectric 
■ functions corresponding to our samples, and for that typically used in the precise prediction of the 

force. The force for our films was found to be 5-14% smaller at a distance of 100 nm between the 
plates. Noise in the optical data is responsible for the force variation within 1%. It is concluded 
^ . that prediction of the Casimir force between metals with a precision better than 10% must be based 

^ I on the material optical response measured from visible to mid-infrared range. 
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I. INTRODUCTION 



The force between electrically neutral metallic plates separated by a small vacuum gap (of micron/submicron di- 
mensions), as predicted by the eminent Dutch physicist H.B. Casimir in 1948, still attracts considerable interest. On 
one hand this interest is curiosity driven since the force is connected with the zero-point fluctuations in vacuum, 
and on the other hand the interest stems from practical applications because modern microtechnology approaches 
the limit where the force starts to influence the performance of microdevices. During the past decade the Casimir 
force was measured with increasing precision in a number of experiments using different techniques such as torsion 
pendulum 2 ., atomic force microscope (AFM) 3 ? 4 -, microelectromechanical systems (MEMS)£&££ and different geomet- 
rical configurations: sphere-plate 2 -^, plate-plate^, crossed cylinders 1 ^. In most cases the bodies were covered with 
gold evaporated or sputter deposited to a thickness of 100-200 nm. 

Relatively low precision, 15%, in the force measurement was reached for the plate-plate configuration? because 
of the parallelism problem. In the torsion pendulum experiment^ the force was measured with an accuracy of 5%. 
In the experiments^^ errors were claimed on the level of 1%. In the most precise up to date experimental^ the 
experimental errors claimed to be as low as 0.5%. 

Comparison between theory and experiment answers an important question: how accurately do we understand the 
\ origin of the force? To make a precise evaluation of the force taking into account real conditions of the experiments 
is equally difficult as to make a precise measurement. In its original form, the Casimir forced given by 



F c (a) = j (1) 

V 7 240 a 4 y ' 

was calculated between the ideal metals at zero temperature. It depends only on the fundamental constants and the 
distance between the plates a. The force between real materials was derived for the first time by Lifshit a 11 ' 12 ' 13 . The 
material properties enter the Lifshitz formula via the frequency dependent dielectric function e(cu). This formula 
became the basis for all precise calculations of the force. Corrections to Eq. ([T]) can be very large especially at small 
separations (< 100 nm) between bodies. The Lifshitz formula accounts for real optical properties of the materials, and 
for finite temperature effects. An additional source of corrections to the force is the surface roughness of interacting 
platesii^. 

In all the experiments mentioned above the bodies were covered with metallic films but the optical properties of 
these films have never been measured. It is commonly accepted^i&ii that these properties can be taken from the 
handbooks tabulated dat a 18 i 19 together with the Drude parameters, which are necessary to extrapolate the data to 
low frequencies. This might still be a possible way to estimate the force, but it is unacceptable for calculations with 
controlled precision . Lamoreaus* 2 ^ was the first who recognized this problem. The reason is very simpl o 20 ' 21 i 22 ' 23 : 
optical properties of deposited films depend on the method of preparation, and can differ substantially from sample 
to sample. 
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Recently analysis of existing optical data for Au was undertaken^! to explore how significant is the effect of variation 
of the optical properties on the Casimir force. It was demonstrated that different sets of the data deviate considerably. 
This variation influences the Casimir force on the level of 5% in the distance range of the most precise experiments. 
Significant sample dependence of the force raises doubts on reported agreement between theory and experiment within 
1% precision^. This is an important issue that has to be further thoroughly investigated. 

In this paper we present the optical properties of Au films of different thickness deposited in two different evaporators 
on silicon and/or mica substrates, unannealed or annealed after the deposition. Moreover, we will discuss the influence 
of measured optical properties of gold films on the precise evaluation of the Casimir force. For the first time the 
characterization of the films was performed over a wide frequency range. It was done ellipsometrically using the 
infrared variable angle spectroscopic cllipsometer in the wavelength range 1.9 - 33 um, and the vacuum ultraviolet 
ellipsometer in the range 0.14 - 1.7 /im. In addition, the film roughness was characterized with atomic force microscopy 
(AFM), and the electron density in the films was estimated from X-ray reflectivity measurements. 

Careful analysis of the data is performed to extract the values of the Drude parameters. It includes joint fits 
of the real and imaginary parts of the dielectric function, or refractive index and extinction coefficient in the low 
frequency range; the Kramers-Kronig consistency of the dielectric function or complex refractive index performed at 
all frequencies. The most important conclusions that follow from this analysis is that the films deposited at the same 
conditions, but having different thicknesses, have considerably different dielectric functions; annealing or change of the 
deposition method showed also influence on the optical properties. At any rate this difference cannot be ignored in a 
precise calculation of the Casimir force. We demonstrate that the optical data typically used for the force evaluation 
in former studies are far away from that found in our samples. The main reason for this deviation is the use of the 
Drude parameters, which correspond to a perfect gold single crystal rather than real polycrystalline films, containing 
a number of different defects. 

The paper is organized as follows. In Sec. [II] we describe preparation and characterizations of Au films and make 
comparison with results known from literature. Analysis of the optical data is presented in Sec. IIII1 where the Drude 
parameters, and uncertainties in these parameters are determined. Calculation of the Casimir force for our samples 
is given in Sec. IIVI Our conclusions are presented in the last Section. 

II. EXPERIMENTAL 

Five gold films were prepared by Au deposition on cleaned (100) Si substrates and freshly cleaved mica. The native 
oxide on Si substrates was not removed; the root-mean-square (rms) roughness was 0.3 nm for the Si substrates, 
while the mica substrate was atomically flat. The first three samples (numbered as 1, 2, and 3) were prepared on Si 
covered first with 10 nm adhesive sublayer of titanium followed by deposition of 400, 200, and 100 nm of Au from 
the source of 99.999% purity. The electron-beam evaporator was used for deposition at a base pressure of 10~ 6 mbar. 
The deposition rate was 0.6 nm/ s. The temperature of the samples was not controlled in the evaporator and it was 
approximately at room temperature. The other two samples were prepared in a thermal evaporator at the same base 
pressure and deposition rate. One film was deposited to a thickness of 120 nm on Si with chromium sublayer (sample 
4). The other film of the same thickness was deposited on mica, annealed at 375 °C (2 hours) and slowly cooled down 
in a period of 6 hours resulting in an atomically flat film (sample 5). 

The Atomic Force Microscope (Veeco Dimension 3100) was used to determine the surface morphology. The rough- 
ness scans are shown in Fig. Q] for all 5 samples. The corresponding rms roughness, w, and correlation lengt h 25 i 26 i 27 £, 
shown in each panel, were obtained as the average values found from multiple scans. The correlation lengths for the 
first 3 samples corresponding to the lateral feature sizes were reported before^. It should be noted that for sample 
4 (120 nm Au/Si) the correlation length is larger than those for the other three films on Si. It can result from 
differences in the evaporation process or due to different adhesive layer. The annealed film on mica had very smooth 
hills and the largest correlation length. 

Optical characterization of the films was performed by J. A. Woollam Co., Inc^. The vacuum ultraviolet variable 
angle spectroscopic ellipsometer (VUV-VASE) was used in the spectral range from 137 nm to 1698 nm at two angles 
of incidence 65° and 75° (±0.01°). The steps in the wavelength A was increased quadratically with A from 1.5 nm 
to 200 nm. In the spectral range from 1.9 am to 32.8 um the infrared variable angle spectroscopic ellipsometer 
(IR-VASE) was used at the same incidence angles in steps <~ A 2 ranging from 1.4 nm to 411 nm. 

From ellipsometry the ratio of p-polarized and s-polarized complex Fresnel reflection coefficients is obtaine d 30 i 31 : 

p= r JL = tan*e lA , (2) 
r s 

where r PiS are the corresponding reflection coefficients, and the angles ^ and A are the raw data collected in a 
measurement as functions of A. All our films can be considered as completely opaque, and can be described by the 
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FIG. 1: (Color online) Surface scans of all films with AFM using the same color scale, 
thickness, rms roughness w, and correlation length £ are shown in each panel (in nm). 
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(3) 



where t? is the angle of incidence and < e >=< e(A) > is the "pseudo" dielectric function of the films. The term 
"pseudo" is used here since the films may not be completely isotropic or uniform; they are rough, and may contain 
absorbed layers of different origin because they have been exposed to air. If this is the case then the dielectric function 
extracted from the raw data will be influenced by these factors. For our films we expect high level of isotropy but 
they can be not uniform in depth. It means that the dielectric function we extract from the data will be averaged 
over the distance of the order of the penetration depth. The roughness and absorbed layer can have some significance 
in the visible and ultraviolet ranges but not in the infrared, where the absorption on free electrons of metals is very 
large. Moreover, the effect of roughness is expected to be small since for all films w is much smaller than the smallest 
wavelength 137 nm. Because the infrared domain is the most important for the Casimir force we will consider < e(A) > 
extracted from the raw data as a good approximation for the dielectric function of the top layer of a gold film. 
The dielectric function is connected with the ellipsometric parameter p for an isotropic and uniform solid as 



e = sin 2 d 



tan 2 ■& 



1 



(4) 



Instead of using e a material is often characterized by the complex refractive index n = n + ik = ^/e, where n is 
the refractive index and k is the extinction coefficient. Both descriptions are equivalent but the noise in the data is 
weighted differently, and it can influence to some degree the values of the Drude parameters (see next Section) . 

Figgure (Ha) shows the experimental results for e"{oS) found via Eq. Q for 3 different films. One can see that 
the IR-VASE gives noisy signal at both ends of the spectral interval. The noise is significant for A > 20 nm but the 
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FIG. 2: (a) Measured e" as a function of frequency u>. (b) The same for |e'|. 



Sample 


A = 1 [im 


A = 5 \xm 


A = 10 jim 


1, 400 nm/Si 

2, 200 nm/Si 

3, 100 nm/Si 

4, 120 nm/Si 

5, 120 nm/mica 


-29.7 + i2.1 
-31.9+i2.3 
-39.1 +i2.9 
-43.8 + i2.6 
-40.7 + il.7 


-805.9 + il85.4 
-855.9 + il95.8 
-1025.2 +i264.8 
-1166.9 +i213.9 
-1120.2 + il78.1 


-2605.1 + il096.3 
-2778.6 + il212.0 
-3349.0 + il574.8 
-3957.2 + il500.1 
-4085.4 + «1440.3 



TABLE I: Dielectric function for different samples at fixed wavelengths A = 1, 5, 10 fim. 

number of points in this range is not large, and the weight of these points for the extraction of the Drude parameters 
or Kramers-Kronig analysis is small. Around the interband transition (minimum of the the curves) the smallest 
absorption is observed for the annealed sample on mica indicating the smallest number of defects in the sample^. On 
the contrary, this sample shows the largest \e'\ in the infrared as one can see in Fig. [TJb). An important conclusion 
that can be drawn from our measurements is the sample dependence of the dielectric function. The log-log scale is 
not very convenient for having an impression of this dependence. We present in Table |T] the values of e for all five 
samples at chosen wavelengths A = 1,5, 10 fim. One can see that the real part of e varies very significantly from 
sample to sample. 

The sample dependence of the dielectric function can be partly attributed to different volume of voids in the 
films as was proposed by Aspnes et al£l. To check this assumption we did standard X-ray reflectivity (XRR) 
measurements22i2i25 f or the 100, 200, and 400 nm Au films on Si, which were deposited at identical conditions. From 
this kind of measurements one can draw information on the density of thin films. For this purpose the Phillips Xpert 
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FIG. 3: X-ray reflectivity (counts) vs the angle of incidence (degrees). All 3 films were deposited at similar conditions. 



diffractometer on the Cu K a radiation line A = 1.54 A was used. The angle between the source and the surface was 
increased from 0.06 to 2 degrees. For hard X-rays the refractive index can be present as n w 1 — S, where S -C 1 and 
we neglect the imaginary part of n since absorption is small. For gold S is given 

S = JL N e " A 2 « 2.70 x lO" 6 -^^ ( -X) , (5) 
2tt 47re mc 2 1 cm" 3 \l A J 

where Z is the atomic number, N a is the number of atoms per unit volume, eo is the permittivity of vacuum, and 
N e = ZN a is the density of electrons. 

Since X-rays refract away from the normal on the surface (refractive index n < 1), there exists a critical angle. 
Below this angle the total reflection occurs. The critical angle 6 C can be related to S as 8 C « \25. Therefore, by 
measuring 6 C one can provide knowledge of the electron density N e . The XRR results are shown in Fig. [3J For very 
small angles the reflectivity decreases (not shown here), which may be due to the beam falling off the sample. Below 
the critical angle the material reflectivity is generally not very well understood, but it is of no concern to us since we 
are only interested in the region around 9 C . The transition region is clearly visible on the graph. Above the critical 
angle the signal drops very fast. The slope depends on the surface roughness. It is clear from the graph that the 
100 nm film (sample 3) has the largest critical angle and, therefore, the largest electron density N e . For bulk gold 
we have N e w 4.67 x 10 24 cm -3 . For our films we found from the critical angles N e « (4.5 ± 0.8) x 10 24 cm" 3 for 
sample 3 and N e ~ (3.6 ± 0.8) x 10 24 cm~ 3 for the two other films. The errors are rather large because N e ~ d\ and 
the curve is not sharp at 9 C . From this measurement we cannot extract quantitative information, but qualitatively it 
agrees with the suggestion of different volume of voids in the films. 



A. Comparison with the existing data 

In the interband absorption region (ui > 2.45 eV) there is significant amount of data obtained by combined 
reflectance and transmittance, ellipsometric spectroscopies on unannealed or annealed thin films or bulk samples 
measured in air or ultrahigh vacuum. For comparison we have chosen the films by Theye 3 - evaporated in ultrahigh 
vacuum (lO" 10 - lO" 11 Torr) in the thickness range 10-25 nm and being well annealed. These films represent the 
bulk-like material. The data were collected by measuring reflectance and transmittance of the films. The Theye data 
became part of the handbook tabled in the range 1 < u> < 6 eV. The second choice is the data by Johnson and 
Christy 37 , These films 25 - 50 nm thick and they were thermally evaporated in vacuum at 10 6 Torr. The data 
were collected for unannealed films from reflection and transmission measurements. As the third choice we took the 
data by Wang at al. 38 . The films were thermally evaporated at pressure 10 -5 Torr. The data were collected with 
spectroscopic ellipsometry for unannealed films of thickness 150 nm. The preparation conditions and the method of 
measurement are similar to that for our films. 
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FIG. 4: Imaginary part of the dielectric function in the interband region obtained in different experiments. The thin solid line 
represents the data by Tlieye 3 -^, the solid circles are the data by Johnson and Christy 37 , and the open circles are the data by 
Wang et al.— . The data of this study are presented by thick solid (sample 3) and thick dashed (sample 2) lines. 

The imaginary part of the dielectric functions, e"(ui), for the chosen experiments are shown in Fig. [4] together with 
our data in the interband range. One can see that in this range of frequencies our data are rather typical. Only the 
Theye films, obtained in ultrahigh vacuum and being well annealed, have significantly larger interband absorption. 
It should be mentioned that the interband absorption correlates with the film thickness: the thicker the film the less 
it absorbs. The smallest absorption is observed for our 200 nm film. Wang et al. 38 deposited their films on bottom 
of Dove prisms and measured the optical response on gold-air and gold-glass interfaces. They have found that the 
interband absorption on the gold-glass interface is larger than that on the gold-air interface. In the former case the 
absorption is close to that observed by Theye. It can be attributed to more dense parking of Au atoms nearby the 
substrate. 

In the infrared and especially in mid- and far-infrared the experimental data are sparse. In this range of frequencies 
both e 1 and e" were measured only in a few studie o 39 ' 40 ! 41 . Padalka and Shklarevskii 41 thermally evaporated the 
films on glass at pressure 10 -5 Torr. The films were not annealed; the thickness of the films was not reported. They 
measured the optical constants in the range 1 < A < 12 [im. Motulevich and Shubin 40 evaporated gold on polished 
glass at pressure ~ 10~ 6 torr. The investigated films were 0.5 — 1 /im thick. The samples were annealed in the 
same vacuum at 400° C . The optical constants n and k were measured by the polarization methods in the spectral 
range 1 — 12 [im. Dold and Mecke did not describe the sample preparation carefully. It was only reported 3 ^ that the 
films were evaporated onto a polished glass substrate, and measured in air by using an ellipsometric technique in the 
range 1.25 — 14 fim. Presumably they were not annealed. These data are included in the handbook 18 , table in the 
corresponding spectral range. 

Fig. [5] shows all the literature low-frequency data and three of our films. The main conclusion that can be drawn 
from this figure is that our films are typical in the sense of optical properties. This is because all the films were 
deposited at similar conditions. The annealed films by Motulevich and Shubin 4 ^ and our annealed sample 5 show the 
largest —e'. 



The Casimir force given by the Lifshitz formula depends on the dielectric function at imaginary frequencies: s(iQ. 
This function cannot be measured directly in any experiment but with the help of the Kramers-Kronig relation it can 
be expressed via the observable function e"(u): 



III. ANALYSIS OF THE DATA 
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FIG. 5: Comparison of existing low-frequency data for gold. Panel (a) gives e" and panel (b) shows — e' as functions of the 
wavelength. The literature data are marked as PM 4 ^, MS^, and DM^ (included in the handbook^). Our data are presented 
for samples 2, 3, and 5. 



The experimental data available for e"(w) are always restricted from low and high frequencies. The low frequency 
cutoff u! cut is especially important in the case of metals. This is because for metals e" is large at low frequencies 
which contribute significantly to e{iC^^. Therefore, an important step in the evaluation of e(iC) is extrapolation of 
the dielectric function e"{u>) to low frequencies lj < Lo cu t, where the experimental data are not accessible. 
At low frequencies the dielectric function of metals can be described by the Drude function 

e{u) = l- I (7) 

which is defined by two parameters, the plasma frequency lu p and the relaxation frequency uj t . Lambrecht and 
Reynaud 42 fixed the plasma frequency using the relation 

where N is the number of conduction electrons per unit volume, e is the charge and m* is the effective mass of 
electron. The plasma frequency was evaluated assuming that each atom gives one conduction electron and that the 
effective mass coincides with the mass of free electron. The bulk density of Au was used to estimate N. The value 
of Up = 9.0 eV found in this way was largely adopted by the community^iSiL2i2iifl. A relatively close value of uj p 
was found by Bennett and Bennett 43 for carefully prepared films deposited in ultrahigh vacuum 4 ^. However, it was 
stressed by these authors that the reflectance for their films was always higher than the values reported for films 
under standard vacuum conditions. The value lj p = 9.0 eV is close to the plasma frequency in a perfect single crystal 
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but the films used for measurement of the Casimir force can contain defects, which are responsible for the reduction 
of Up. The most important defects are small voids, which were observed in gold films with transmission electron 

45 46 

microscop y ' . 

Special investigation of the influence of the defects was undertaken by Aspnes et al^, where it was stressed 
that films grown at different conditions have considerably different e" above the interband transition. The spectra 
qualitatively differ only by scaling factors as one can see in Fig. [4j The scaling behavior of the spectra was attributed 
to different volumes of voids in the films prepared by different methods 32 . Different kind of defects will also contribute 
to the scattering of free electrons changing the relaxation frequency uj t . Many researches stressed that the conduction 
electrons are much more sensitive to slight changes in the material structur o 32 i 36 ' 44 i 47 . It is well known, for example, 
that the resistivity of a film can be significantly larger than the resistivity of the bulk material. At any rate, so far we 
have to conclude that both Drude parameters must be extracted from the optical data of the films, which are used 
for the Casimir force measurement. Below we describe a few ways to extract these parameters from our data. 



A. Joint fit of e' and e' 



Separating the real and imaginary parts of the Drude function ([7]) one finds for s' and e": 



£ » = 1- 



1 i ? ' 



W (tjJ 2 + W 2 ) 



(9) 



These equations have to be true below the interband transition Tiuj < 2 Ah eV (A > 0.5 /im) 36 , but because this 
transition is not sharp one has to do analysis at lower frequencies. Practically Eqs. ((9]) can be applied at wavelengths 
A > 2 /im that coincides with the range of the infrared ellipsometer. In this range Eqs. © can be compared with the 
optical data for both e'(u>) and Minimizing deviations between the data and theoretical expectations one can 

find the Drude parameters io p and ui T . For example, for sample 3 (100 nm, Au/ Si) we found uj p — 7.79± 0.01 eV and 
u) T = 48.8 ± 0.2 meV . Similar calculations for annealed sample 5 (120 nm, Au/mica) gave oj p = 8.37 ± 0.03 eV and 
u) T = 37.1 ± 0.5 meV. The statistical uncertainty of the Drude parameters was found using a \ 2 criterion for joint 
estimation of two parameters 4 ^. For a given parameter the error corresponds to the change Ax 2 = 1 when the other 

parameter is kept constant. For all samples the values of the parameters are collected in Table [II] 

We found the parameters also in a slightly different way. One can fit the complex refractive index h(oj) = ^/e(u) 
instead of the dielectric function. In the Drude range, where nearly all absorption is due to free electrons in the metal, 
n(w) and k(ui) behave as 



n(uj) 



V2oj t y/l + y 2 



1 -1/2 



fc(w) 



V2Wt y/T 



1/2 



(10) 



where y — lo/lu t . Then we can minimize deviations for n(u>) and k(u>). The corresponding parameters for sample 3 
are uj p = 7.94 ± 0.01 eV and uj t = 52.0 ± 0.2 meV. For the annealed sample 5 they are uj p = 8.41 ± 0.02 eV and 
u) T = 37.7 ± 0.4 meV. It can be noted that within the statistical errors the parameters for the annealed film are the 
same as those found by joint fit of e' and e" . However, for sample 3 this is not the case. 

The Drude parameters should be the same in both cases but some difference can appear due to the smaller con- 
tributing weight of low frequencies when we perform minimization for n and k than that from the minimization of 
e' and e". Figure [6] shows the data for e'(X) and e"(X) (solid lines) and the best Drude fit (dashed lines) found by 
minimization of deviations for e' and e" . Panel (a) corresponds to sample 5 and panel (b) shows the data for sample 3. 
The film on mica is described well by the Drude dielectric function, but the film on Si demonstrates some deviations 
from the Drude behavior at A < 15 fim. It looks like an additional absorption band. In principle, the anomalous skin 
effect can be responsible for absorption in this range, but we found that this effect can be observable only at smaller 
lu t . Additional absorption in the Drude range is often observed due to the tail of the interband transition, but this is 
hardly the case because the wavelength is too large. 

In absence of information on the origin of this absorption band we did phenomenological analysis of our data with 
the dielectric function, which includes an additional Lorentz oscillator: 



e (w) = e D (uj) + 



Swl 



(11) 



where Sd is the Drude function ([7]), S is the dimensionless oscillator strength, lu and 7 are the central frequency and 
width of the band, respectively. Fig. [7] shows the difference Ae" = e" — e" D as a function of the wavelength A. The 
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FIG. 7: Additional absorption band for samples 2, 3, and 5 as a function of the wavelength. 



smallest absorption is realized for the annealed film on mica. It becomes larger for the unannealed 100 nm film on 
Si, and increases further for the 200 nm film on Si. On the other hand, the central wavelength of the band and its 
width do not change significantly from sample to sample showing the common origin of the band for different samples. 
This situation can be expected if the samples differ only by the density of defects of the same kind. However, without 
additional experimental information we cannot specify the exact nature of these defects. 
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FIG. 8: Dielectric function e' as a function of uj for sample 3. Thick line is the experimental curve extrapolated to the frequencies 
lower than 0.038 eV according to the Drude model. Thin line is the same function calculated from the K-K relation (|12p . Inset 
shows the behavior of e' inside of the dotted rectangular in the linear scale. 

B. Determination of the parameters using the Kramers-Kronig relation 

One of the disadvantages of the ellipsometric determination of the dielectric function is that the method does not 
maintain Kramers-Kronig (K-K) consistency. Therefore, it is important to check the K-K relations for our data. To 
do the K-K analysis we also have to extrapolate the dielectric function outside of the measured frequency range. For 
metals extrapolation to low frequencies is based again on the Drude dielectric function ([7]). It means that together 
with K-K consistency we will find the Drude parameters for which this consistency is the best. This method for 
determining the Drude parameters was used recently^ for analysis of the optical properties of gold samples, which 
are available in the literature. 

The K-K relation expresses e'(uj) as integral over all frequencies of e"(uj) : 

2 7 , xe" (x) , . 

e' (w - l = -P dx^^, (12) 

7T J X — UJ 



where P means the principal part of the integral. To use this relation we have to define s"(uj) at all frequencies. Below 
the low- frequency cutoff we define it using Eq. © . Above the high frequency cutoff (9 eV for our data) we enlarge the 
frequency range by using the handbook datai& in the range 9 < ui < 100 eV and above 100 eV s"(uj) is extrapolated 
as A/uj 3 . The constant A is determined by matching e" at uj = 100 eV. The principal part of the integral in (fT2"|) 
is calculated in the same way as ir*2i. e'(uj) was calculated in the frequency range 0.01 < uj < 100 eV as a function 
of the Drude parameters. This function was compared with the measured one and minimization of deviations gave 
us the values of the parameters. The points outside of the measured frequency range (uj < 0.038 eV) were compared 
with the prediction based on Eq. (J9j> . 

For sample 3 it was found uj p = 7.80 eV and uj t = 47.9 meV . The experimental function e' (uj) continued to lower 
frequency according to the Drude model and its prediction based on the relation (fT2"]l are shown in Fig. [5] by thick 
and thin lines, respectively. The agreement between the curves is rather good. At high frequencies where |e'| ~ 1 the 
logarithmic scale is not convenient for comparison. Instead we present in the inset the region bounded by the dotted 
rectangular in the linear scale. In this range the positions of the peaks are reproduced very well, but their magnitude 
is slightly different. This is because we used the data for e" above 9 eV from the handbook, which did not match 
precisely to those for our sample. Similar situation is realized for all the other films. The Drude parameters for all 
samples are presented in Table HU 

Alternatively one can use the experimental extinction coefficient k(uj) to get the refractive index n(uj) using the 
K-K relation between n and k: 

oo 

n(uj)-l = -P fdx^H,. (13) 

7T J X z — UJ Z 
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FIG. 9: Refractive index n as a function of ui for sample 3. Thick line is the experimental curve extrapolated to u) < 0.038 eV 
according to the Drude model. Thin line is the same function calculated from the K-K relation (|13[) . The inset shows the 
behavior of n in the linear scale in the range 0.5 < uj < 50 eV. 



At frequencies for which the experimental data are not accessible we define k(u>) — Im^/ e{uS), where s(lo) continued 
to low frequencies according to (|7|) and to high frequencies as e(u>) = 1 — ^ju) 2 + iA/ui 3 . The constant A again is 
chosen by matching e" at the highest frequency u) — 100 eV. 

In contrast with the K-K relation (|12p now we cannot present the contribution of low frequencies to the dispersion 
integral (fT3|) in the analytic form. The Drude representation (fT0|) for k(u>) at uj -C uj p is 



k(w) = 



\/2u> r 



1 



1 



1/2 



(14) 



where y = lu/lo t . This function of y was approximated by the forth order polynomial in y in the range < y < 2. With 
this polynomial the contribution of low frequencies to (| 13[) was found analytically. Minimization of deviations between 
the experimental values of n(ui) and the theoretical predictions via (|13p gave the values of the Drude parameters. 

For sample 3 it was found ui p = 7.84 eV and u> T = 47.4 meV. The experimental and predicted curves for n(ui) are 
shown in Fig. [9] by thick and thin lines, respectively. The inset shows the same functions in the linear plot in the 
range 0.5 < u> < 50 eV. Again we have a reasonable agreement of the experiment and prediction on the basis of the 
K-K relation. 



The values of the Drude parameters found by different methods described above are collected in Table ITT] The 
statistical errors in the parameters are rather small. They do not depend significantly on the method and vary only 
slightly from sample to sample. These errors are 0.01 — 0.03 eV for lo p and 0.2 — 0.5 meV for u> T . As one can see 
from Table |TT] the values found by different methods do not agree with each other within the statistical errors. This is 
because each method treats noise in the data differently. We cannot give preference to any specific method. Instead, 
we average the values of the parameters determined by different methods, and define the root-mean-square (rms) 
error of this averaging as uncertainty in the parameter value. The averaged parameters and rms errors are given in 
the last column of Table [Til Samples 1 and 2 have similar Drude parameters, which cannot be resolved within the 
discussed uncertainty, but all the other samples arc clearly different. 



IV. SENSITIVITY OF THE CASIMIR FORCE TO THE OPTICAL PROPERTIES OF GOLD FILMS 



In this Section we are going to discuss the Casimir force without temperature or roughness corrections in order to 
concentrate on the influence of the material optical properties on the force. In this case the force between two similar 



12 



Q Q TTI T"\l p 


Pa T" a vnpi-pT- 
J. aiciiiicici 


Tnint r' f" 

U *J ill Li c d c 


Inint t? 1c 


K-K e' 


K-K n 


A irpv^* CCP 


1 


^ p [eV] 


6.70 


6.87 


6.88 


6.83 


6.82±0.08 


400 nm/Si 


uj T [meV] 


38.4 


43.3 


40.2 


39.9 


40.5±2.1 


2 


Ulp 


6.78 


7.04 


6.69 


6.80 


6.83±0.15 


200 nm/Si 


LU T 


40.7 


45.3 


36.1 


36.0 


39.5±4.4 


3 


Wp 


7.79 


7.94 


7.80 


7.84 


7.84±0.07 


1UU 71771, j o% 


UJ T 


4o.o 


oz.u 


A1 O 


A1 A 


ylQ 0-1-0 1 


4 


Ulp 


7.90 


8.24 


7.95 


7.90 


8.00±0.16 


120 nm/Si 


LU T 


37.1 


41.4 


35.2 


29.2 


35.7±5.1 


5 


Up 


8.37 


8.41 


8.27 


8.46 


8.38±0.08 


120 nm/mica 


U-r 


37.1 


37.7 


34.5 


39.1 


37.1±1.9 



TABLE II: The Drude parameters determined by different methods described in the text. In all cases the statistical errors in 
the parameters are on the same level: 0.01 — 0.03 eV for lu p and 0.2 — 0.5 meV for U) T . The last column shows the values of the 
parameters averaged on different methods and the corresponding rms errors. 



parallel plates can be calculated using the Lifshitz formula 13 : 

co oo 

Fpp (a) = -gla / d W d 11 ko E i.^e^afco ' ( 15 ) 

where q is the wave vector along the plates (q — |q|). The formula includes the reflection coefficients for two 
polarization states [i = s and fi = p which are defined as 

_ k - k x _ e «) fc - fci 

rs ~fc + fei' rp " £ (*C)fco + fci 1 j 

with fco and k\ being the normal components of the wave vector in vacuum and metal, respectively: 

fco = VCV^ + g 2 , fci = V^0(W+7- (17) 



A. Dielectric function at imaginary frequencies 



To evaluate the force with the Lifshitz formula one has to know the dielectric function of the material at imaginary 
frequencies e(i(), which is calculated via e"(w) according to Eq. (JSJ). Therefore, first we have to calculate e(iQ using 
our optical data. For this purpose let us present e(iC) as 



£«) = 1 + e cut (i() + e exp (i(), 



(18) 



where e cu t is calculated with the extrapolated s"{u>) in the unaccessible frequency range lo < uj cu t, while e exp is 
calculated using the experimental data for e"(w), according to the formulas: 



Zcut{K) 



Loe"{u) 
' Q 2 +uj 2 



rp(*C) 



dio 



' C 2 + OJ 2 



(19) 



Strictly speaking, we included in e exp the extrapolation to high frequencies, uj > 100 eV, but this is justified because 
these frequencies do not play very significant role. The high- frequency extrapolation of e"(u>) was done in the same 
way as in Sec. lIIIBl so as the calculation of the integral for e exp . For the low-frequency extrapolation the Drude model 
([9]) was used. In this case the integral for e CMt can be found analytically, and it yields 



ir( 2 



tan 



--tan 



(20) 



Note that there is no singularity here at £ = u> T . 

It was already stressed that for metals e cut gives an important contribution to the dielectric function. Of course, it 
depends on the value of uj cu t- For all previous data this value was around 0.125 eV. In this study u> C ut = 0.038 eV is 
about 3 times smaller, but still the contribution of e cut is significant. It can be seen from Fig. 1101 where the relative 
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FIG. 10: Relative contribution of low (u < uj cu t) and high (w > u> C ut) frequencies to the dielectric function e(i£) for our sample 
3 (thin lines) and for the handbook data 18 (thick lines). 

values £ C ut(iC) I ' e (*C) ancl £ezp(*C)/ £ (*C) are presented for sample 3 as thin lines. For comparison in the same plot we 
showed (thick lines) e cu t and e exp calculated with the handbook data, and extrapolated to lo < 0.125 eV with the 
Drude parameters io p = 9.0 eV and uj t = 35 meV. It should be stressed that for our film the contribution from the 
extrapolated region, e cu t, dominates at £ < 0.2 eV, while for the handbook data it dominates up to £ = 4 eV. This 
is the result of reduced Lu cut for our data. It means that the calculations based on our data are more reliable because 
a smaller part of e(££) depends on the extrapolation. 

The total dielectric functions £j(i£) {i 1S the number of the sample) are presented in Fig. [TIT a) for samples 2, 3, 
and 5. The results for samples 1 and 4 are not shown for clarity. The thick solid line represents £o(«£), which is 
typically used for the Casimir force calculation. In this case the integral for e exp is evaluated for e"(uj) taken from the 
handbook™, where the cutoff frequency is u> cu t = 0.125 eV. The contribution of low frequencies, £ CM t, is calculated 
with the Drude parameters u p = 9.0 eV and oj t = 35 meV—. In what follows we are using £o(?£) as a reference case. 

Important contribution to the Casimir force comes from the imaginary frequencies around the characteristic fre- 
quency Cc/i = c/2a. This frequency is in the range 0.1 < £ c /j < 10 eV when the distance is in the most interesting 
interval 10 nm < a < 1 \xm. The frequency range in Fig. Illf a). (b) was chosen accordingly. The logarithmic scale 
in Fig. [TlTa') does not give the feeling of actual difference between the curves. The relative change in the dielectric 
function [eo(i£) — £i(i£)] / £ o(*£) demonstrates much better the significance of actual optical properties of the films. 
This change is shown in Fig. [TiTb) for all our films. The deviation of £i{iQ from the imaginary material with the 
dielectric function £o (*C) > which is described by the bulk Drude parameters and the handbook optical data, is signifi- 
cant at all important frequencies and for all samples. Of course, this deviation will be translated into the change in 
the Casimir force. The curves in Fig. [TlT b) were calculated for the middle values of the Drude parameters in the last 
column of Table ITT1 The uncertainty of these parameters is responsible for uncertainty of £i(iQ. It is especially large 
for samples 2 and 4. We will discuss it later in connection with the uncertainty of the force. 



To calculate the Casimir force the dielectric functions for all samples, £j(z£), were found numerically in the frequency 
range 0.01 eV < ( < 100 eV. At lower frequencies, £ < 0.01 eV, they were extrapolated according to the Drude 
model. At higher frequencies, £ > 100 eV, we extrapolated with the function £j(«C) = 1 + A-i/C' ''■> where for each 
sample the constant Ai was chosen to match the value of £i(i() at £ = 100 eV. 

It is convenient to calculate not the force itself but so called the reduction factor rj pp , which is defined as the ratio 
of the force to the Casimir force between ideal metals: 



B. The Casimir force 



(21) 
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FIG. 11: (a) The dielectric function as a function of the imaginary frequency (. The thin dotted, dashed, and solid lines 
correspond to samples 2, 3, and 5, respectively. The thick line gives eo(iC) that is typically used for the Casimir force 
calculations, (b) Relative deviation of the dielectric function of i-th sample, from eq. 



p- 0.5 




FIG. 12: Reduction factor r] as a function of the separation a for samples 1,3, and 5. The thick line shows the reference result 
calculated with eo(iQ- 



For convenience of the numerical procedure one can make an appropriate change of variables in Eq. (|15p so that the 
reduction factor can be presented in the form 



'w<«') = £ E M-$rp (22) 



2tt 4 



where the reflection coefficients as functions of x and y arc defined as 



1 + s' p e{iC, ch xy) + s ' 

with 



1-s ejiCchxy) - s 

r p = ~T7 r~ — i (23) 



s = Vl + x 2 [e^Cchiy)-!], Cch = ^. (24) 
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FIG. 13: Relative deviation of the reduction factors for different samples from the reference curve r/o(a), which were evaluated 
using the handbook optical datai^ and the Drude parameters uj p — 9 eV, lu t = 0.035 eV. 



The integral ([22]) was calculated numerically with different dielectric functions £i(«C) with a precision of 1CP 6 . The 
results are presented in Fig. [T^] for samples 1, 3, and 5. The reference curve (thick line) calculated with £o(«C) is al so 
shown for comparison. It represents the reduction factor, which is typically used in the precise calculations of the 
Casimir force between gold surfaces. One can see that there is significant difference between this reference curve and 
those that correspond to actual gold films. To see the magnitude of the deviations from the reference curve, we plot 
in Fig. [13] the ratio (r/o — r?i)/??o as a function of distance a for all five samples. 

At small distances the deviations are more sensitive to the value of uj p . At large distances the sample dependence 
becomes weaker and more sensitive to the value of ui T . For samples 1 and 2, which correspond to the 400 nm and 200 
nm films deposited on the Si substrate, the deviations are especially large. They are 12-14% at a < 100 nm and stay 
considerable even for the distances as large as 1 am. Samples 3, 4, and 5 have smaller deviations from the reference 
case but even for these samples the deviations are as large as 5-7%. 

We calculated also how uncertainty in the Drude parameters influences the uncertainty in the force. For this purpose 
we calculated the reduction factor r\ at the borders of the error intervals: [lo p + Alu p , lu t ] and [o> p , lu t — Aw T ], where 
Auj p and Aw T are shown as the errors in the last column of Table HH The results were compared with n calculated 
with the middle values of the parameters [lo p , uj t \. The maximal deviations were found for sample 4. The relative 
deviations for this sample, An/rj, are shown in Fig.[T3]as functions of the separation a. These deviations, An/77, are 
defined as 



An r/(ujp, cj t , a) — n(u> p + Su> p , to T + Slu t , a) 
n n(aj p ,u T ,a) 



(25) 



where the variations of the plasma and relaxation frequencies, 5u> p and 5uj t , give the maximal effect on the reduction 
factor when 8lo p = Auj p and 5lj t = — Auj t . 

In most of the experiments the force was measured between a gold covered sphere and a plate. In the case of the 
sphere-plate interaction we define the reduction factor as 

F sp (a) ir 3 Rhc 

Ma) = W ^ (a) = ~^' (26) 

where R is the sphere radius. The expression convenient for numerical calculations in this case is the following: 

1 00 

vs P (a) = ~ J dx J d yy 2 ln I 1 - r l^ v ) > ( 2? ) 



where the reflection coefficients arc defined in (|23l) . Although in the case of sphere and plate the reduction factors 
are different from those between two plates, the qualitative behavior with the separation is the same. The relative 
deviations (770 — ^O/^o for the sphere and plate are very close to those shown in Fig. 1131 and we do not show this 
additional plot. 
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FIG. 14: Relative variation of the reduction factor (see definition in the text) for the Drude parameters at the sides of the error 
intervals. The results are presented for sample 4. 



In this paper we analyzed the optical properties of gold films deposited with different techniques on silicon or mica 
substrates. The optical responses were measured ellipsometrically in a wide range of wavelengths. In the mid-infrared 
our data cover the range of wavelengths up to 33 /im, which is larger than in all previous studies (~ 10 fim), where 
both the real and imaginary parts of the dielectric function were measured. The data unambiguously demonstrate 
the sample dependence of the dielectric function. It was found that thicker films are inferior in the optical sense 
than thin (but opaque) films. It is probably connected with more dense parking of atoms nearby the substrate. The 
X-ray reflectivity measurements qualitatively support this interpretation. We also observed some difference between 
samples deposited in the electron-beam evaporator or in the thermal evaporator, but it is possible that this difference 
is connected with different sublayers (titanium or chromium, respectively) on Si substrates. 

The Casimir force depends on the dielectric function at imaginary frequencies, e(i£), which is expressed via the 
measurable function e"{uj) by the dispersion relation For metals low frequencies play a very important role in 
this relation. However, in experiments there is a low frequency cutoff, co cu t, so that the data at smaller frequencies are 
not accessible. To make precise evaluation of the force, one has to extrapolate the measured e"(u>) to the frequencies 
u) < oJcut- For noble metals like gold it can be done with the help of the Drude model. The parameters of the Drude 
model Up and uj t have to be extracted from the optical data at lu > Lu cut . In this paper we found the Drude parameters 
for our films using joint fit of e'(u>) and s"(lu) and independently using the Kramers-Kronig consistency of the data. 
As an alternative we did the same analysis for the refraction index and the extinction coefficient. The results collected 
in Table |TT] exhibiting significant variation of the parameters from sample to sample. Moreover, the parameters are 
slightly different when different methods to fit the data are used. This is because the noise in the data are weighted 
differently for different methods of parameter determination. We defined the values and the errors in the parameters 
as the averaged values over different methods and the corresponding rms errors, respectively. 

The contribution of the extrapolated part of s"(u>) to is considerably smaller for our data than for the 

handbook data as one can see in Fig. [TU1 The reason is that in our case the cutoff frequency uj cut = 0.0378 eV is 
smaller than that for the handbook data, where u> cu t = 0.125 eV. Nevertheless, this contribution is still significant 
and has to be carefully considered. As a reference curve for e(iQ we have chosen £o(*C)j which was calculated with the 
handbook data for u > 0.125 eV, and extrapolated for smaller frequencies with the Drude parameters uj p = 9.0 eV 
and Lij T = 35 meV. This curve is typically used for the calculation of the Casimir force. For our films we have found 
that £i(iC) was always smaller than the reference curve, and the relative deviation is larger than 15% at £ ~ 1 eV 



Indeed, the Casimir force evaluated for our films is considerably smaller from that calculated with the reference 
curve £o(iC)- The smallest deviation is realized for sample 5 (120 nm Au/mica) and estimated to be 4-5% in the 
distance range a < 200 nm. For thicker films (sample 1 or 2) it can be as large as 14% (see Fig. [13]) . At large 
separations the sample dependence becomes weaker, but it is above 1% even at a = 1 /im. 

Our main conclusion is that actual optical properties of the materials used for the measurement of the Casimir 



V. 



CONCLUSIONS 



(see Fig.QHb)). 
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force are very important for comparison between theory and experiment. These properties have to be measured for 
the same materials, and not taken from a handbook for a material of the same chemical nature but possibly of very 
different microstructure. It is concluded that prediction of the Casimir force with a precision better than 10% must 
be based on the material optical response measured from visible to mid-infrared range. 



* Corresponding author:v. b.svetovoy@ewi.utwente.nl 

1 H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948). 

2 S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997); 81, 5475 (1998). 

3 U. Mohideen and A. Roy, Phys. Rev. Lett. 81, 4549 (1998); A. Roy, C.-Y. Lin, and U. Mohideen, Phys. Rev. D 60, 111101(R) 
(1999). 

4 B. W. Harris, F. Chen, and U. Mohideen, Phys. Rev. A 62, 052109 (2000). 

5 H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, and F. Capasso, Science 291, 1941 (2001); Phys. Rev. Lett. 87, 
211801 (2001). 

6 R. S. Decca, D. Lopez, E. Fischbach, and D. E. Krause, Phys. Rev. Lett. 91, 050402 (2003). 

7 R. S. Decca, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, D. Lopez, and V. M. Mostepanenko, Phys. Rev. D 68, 116003 
(2003). 

8 R. S. Decca , D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepanenko, Ann. Phys. 318 37 
(2005). 

9 G. Bressi, G. Carugno, R. Onofrio, and G. Ruoso, Phys. Rev. Lett. 88, 041804 (2002). 

10 T. Ederth, Phys. Rev. A 62, 062104 (2000). 

11 E. M. Lifshitz, Zh. Eksp. Teor. Fiz. 29, 94 (1956) [Sov. Phys. JETP 2, 73 (1956)]. 

12 I.E. Dzyaloshinskii, E.M. Lifshitz and L.P. Pitaevskii, Advances in Physics 38, 165 (1961). 

13 E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2 (Pergamon Press, Oxford, 1980). 

14 G. L. Klimchitskaya and Yu. V. Pavlov, Int. J. Mod. Phys. A 11, 3723 (1996). 

15 C. Genet, A. Lambrecht, P. Maia Neto and S. Reynaud, Europhys. Lett. 62, 484 (2003). 

16 P. A. Maia Neto, A. Lambrecht and S. Reynaud, Phys. Rev A 72, 012115 (2005). 

17 G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Phys. Rev. A 61, 062107 (2000). 

18 Handbook of Optical Constants of Solids, edited by E. D. Palik (Academic Press, 1995). 

19 J. H. Weaver, C. Krafka, D. W. Lynch, and E. E. Koch, Optical Properties of Metals, Part II, Physics Data No. 18-2 
(Fachinformationszentrum Energie, Physik, Mathematik, Karsruhe, 1981). 

20 S. K. Lamoreaux, Phys. Rev. A 59, R3149 (1999). 

21 S. K. Lamoreaux, Phys. Rev. Lett. 83, 3340 (1999). 

22 V. B. Svetovoy and M. V. Lokhanin, Mod. Phys. Lett. A 15, 1013 (2000). 

23 V. B. Svetovoy and M. V. Lokhanin, Mod. Phys. Lett. A 15, 1437 (2000). 

24 I. Pirozhenko, A. Lambrecht, and V.B. Svetovoy, New J. Phys. bf 8 238 (2006). 

25 P. Meakin Phys. Rep. 235, 1991 (1994). 

26 J. Krim and G. Palasantzas, Int. J. of Mod. Phys. B 9, 599 (1995). 

27 Y. -P. Zhao, G. -C. Wang, and T. -M. Lu, Characterization of amorphous and crystalline rough surfaces-principles and 
applications, (Experimental Methods in the Physical Science Vol. 37, Academic Press, 2001). 
P. van Zwol, G. Palasantzas, and J. Th. M. De Hosson, Appl. Phys. Lett. 91, 144108 (2007). 



http://www.JAWoollam.com 



R. M. A. Azzam and N. M. Bashara, Ellipsometry and Polarized Light (Amsterdam, North Holland, 1987). 



29 
30 

31 H. G. Tompkins and W. A. McGahan, Spectroscopic Ellipsometry and Reflectometry (New York, Wiley, 1999). 

32 D. E. Aspnes, E. Kinsbron, and D. D. Bacon, Phys. Rev. B 21, 3290 (1980). 

33 S. K. Sinha, E. B. Sirota, S. Garoff, and H. B. Stanley, Phys. Rev. B 38 2297 (1988). 

34 A. Braslau, P. S. Pershan, G. Swislow, B. M. Ocko, and J. Als-Nielsen, Phys. Rev. A 38 2457 (1988). 

35 J. Als-Nielsen, Handbook of synchrotron radiation (North Holland New York, 1991). 

36 M.-L. Theye, Phys. Rev. B 2, 3060 (1970). 

37 P. B. Johnson and R. W. Christy, Phys. Rev. B 6, 4370 (1972). 

38 Yu Wang et aL, Thin Solid Films 313, 232 (1998). 

39 B. Dold and R. Mecke, Optik, 22, 435 (1965). 

40 G. P. Motulevich and A. A. Shubin, Zh. Eksp. Teor. Fiz. 47, 840 (1964) [Soviet Phys. JETP 20, 560 (1965)]. 

41 V. G. Padalka and I. N. Shklyarevskii, Opt. Spektroskopiya 11, 527 (1961) [Opt. Spectry. (USSR) 11, 285 (1961)]. 

42 A. Lambrecht and S. Reynaud, Eur. Phys. J. D 8, 309 (2000). 

43 H. E. Bennett and J. M. Bennett, in Optical Properties and Electronic Structure of Metals and Alloys, edited by F. Abeles 
(North- Holland Publ., Amsterdam, 1966). 

44 J. M. Bennett and E. J. Ashley, Appl. Opt. 4, 221 (1965). 

45 J. A. Lloyd and S. Nakahara, J. Appl. Phys. 48, 5092 (1977). 

46 J. A. Lloyd and S. Nakahara, J. Vac. Sci. Technol. 14, 655 (1977). 

47 N. Kaiser, Appl. Opt. 41, 3053 (2002). 



K. Hagiwara et al., Phys. Rev. D 66, 010001 (2002) . 



